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Abstract 



< 

. We present a complete three-loop calculation of the electric dipole moment 

of the u and d quarks in the Standard Model. For the d quark, more relevant 
for the experimentally important neutron electric dipole moment, we find 
1^ I cancellations which lead to an order of magnitude suppression compared with 

•y^ ■ previous estimates. 

^ ■ The electric dipole moment (EDM) of the neutron, dn-, has been subject of experimental 

, searches for almost a half century [||]. Discovery of a non-zero dn would be direct evidence 
^ \ of violation of both time reversal symmetry and parity. If the symmetry under simulta- 

\ neons charge conjugation, space inversion, and time reversal (CPT) is exact, T violation 
is equivalent to CP violation (see also a discussion in 0). Most theories put forward to 



explain CP violation observed in neutral kaon decays, or to address the question of matter- 
^ , antimatter asymmetry in the universe, predict the existence of dn at some level. The current 
experimental upper bound 



^ : rf„ < 1.1 X 10"^^ e cm (1) 

stringently constrains physics both within the standard model (possible non-perturbative 
effects, such as the so-called 6'-term 0) and beyond it, such as supersymmetry and super- 
gravity 0, fourth generation of fermions P], CP violation in the Higgs sector etc. (for 
detailed reviews see |IO| , pd]| ). It also severely constrains the W boson EDM |T^. 



The next generation of experiments is going to search for dn with an accuracy improved 
by several orders of magnitude. For example, a new technique based on storing polarized 
ultra cold neutrons and a polarized gas of ^He in superfiuid ^He is expected to permit a 
measurement of dn with the sensitivity of 

a{dn) = 4 X 10-2^ e cm (2) 



*Address after September 1997: Physics Department, Brookhaven National Laboratory, Upton, 
New York 11973. 
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in a one year run [|T3 

In the electroweak sector of the standard model a complex phase of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix violates CP and, via loop effects, induces an EDM for all 
non-selfconjugate particles with spin. In the case of the neutron, an EDM can be generated 
by electroweak interactions in various ways. Penguin diagrams with a W exchange between 
constituent quarks offer one possibility. Such effects were estimated to give contributions to 
dn which could perhaps be as large as 10"^" e cm WM. This estimate is controversial and 



involves a large uncertainty; see e.g. a discussion in A more conservative estimate of 
the long distance standard model effects is 2 x 10~^^ e cm . 

Another possibility, and the main focus of our work, is dn induced by an EDM of the up 
and down quarks (du^d) inside the neutron. Using SU(6) wave functions one finds 

4 1 

dn = -dd - -du- (3) 

Quark EDM cannot be generated in the standard model at the one loop level because 
the relevant amplitudes do not change the quark flavor and each CKM matrix element is 
accompanied by its complex conjugate; no T-violating complex phase can arise. At the two 
loop level individual diagrams have complex phases (see Fig. |1|) and contribute to the EDM 
|T7HT9|. However, the sum over all quark flavors in the intermediate states leads to the 
vanishing of the EDM at two loops |^ (a similar effect for leptons was found in ref. [pT[]). 



This interesting and seemingly accidental cancellation was also analyzed in |2^j2^ . 

The fact that the quark EDM appears only at the three loop level in the standard model 
greatly complicates theoretical estimates. The largest effect is due to the exchange of two 
W bosons and one gluon; this was flrst discussed by Shabalin in [^] in a model with quarks 
carrying integral charges and 1. In that model the charged quark EDM was found to 
be suppressed by at least 6 inverse powers of the W mass, if one assumes quarks to be 
light. The case of quarks with physical, fractional values of electric charges was examined in 



It was found that the strong cancellation found earlier was an artifact of the assumed 



pattern of charges and quark EDM was found at the level of (quark masses) ^/M^ with 
further enhancements by logs of quark mass ratios; in fact, in this order the EDM of the 
quark is proportional to the electric charge of its isospin partner — so that the EDM of a 
charge 1 quark in Shabalin's model is suppressed. 



The analysis of ref. ||2^, the most complete study to date, was carried out in the frame- 
work of the leading logarithmic approximation within the effective Fermi theory under an 
assumption that all quarks are lighter than the W boson. Since it turns out that the top 
quark outweighs the W boson, this approach is no longer valid. In fact, the large mass of 
the top quark signiflcantly enhances some electroweak processes (e.g. the decays Z bb 
and b — > 57, as well as — mixing). In view of the anticipated improvement of the ex- 
perimental accuracy a detailed evaluation of the standard model contributions to the quark 
EDM is clearly warranted. Here, we report the complete results of such an analysis. 

The quark matrix element of the electromagnetic current, J^™, can be written as (with 
q = p' — p) 
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The form factors at = give the electric charge, anomalous magnetic moment, electric 
dipole moment, and anapole moment in units of e = |e| (e.g. for a down quark -Fi(O) = —1/3), 
so that dq = 6^3(0). 

The three-loop calculation of the light quark EDM is facilitated by the clear hierarchy of 
masses niu^d <^ ^ri^ <C -C mf, -C Mw -C mj. On the other hand, as will be seen shortly, 
the relevant mass ratios are not large enough for their logarithms to be the dominant effect. 
Therefore, in addition to the logarithmic terms we also need the non- logarithmic constants. 
A systematic expansion of Feynman diagrams in terms of mass ratios and their logs is possible 
in the framework of asymptotic expansions (for a recent review see ^^). The problem is 



slightly complicated by the many mass scales present in each diagram. The details of our 
procedures will be described elsewhere. Here we only remark that the calculation technically 
resembles the two-loop electroweak corrections to the muon anomalous magnetic moment 
26| , apart from the presence of the third loop. 

The electric dipole moment in the standard model has an important feature which sim- 
plifies its calculation. Namely, if the source of CP violation is in the complex phase of the 
CKM matrix, the EDM vanishes if any two up-type or down-type quarks have equal masses. 



This circumstance has been taken advantage of in the Fermi theory calculation in ref. ^2 
Under the assumption that all quarks are lighter than the W boson it was possible to replace 
the W propagators by effects of large virtual momenta which could "feel" the exact 

structure of these propagators are, to first approximation, independent of quark masses, 
and therefore their sum gives no contribution to the EDM. The resulting effective theory 
diagrams are shown in Fig. |[ 

The situation is actually more complicated since the top quark is heavier than the W 
boson. It turns out that there are 6 topologies of diagrams with the top quark which 
contribute to the EDM: 3 for the down quark (Fig. |^) and 3 different ones for the up quark 
(Fig. ^). The diagrams in Figs. ^ and ^ have to be understood in the following way: the W 
propagator connected to light quark lines only can be contracted to a point and replaced 
by the other W propagator, connected to the top quark and indicated explicitly in 

the figures, is treated differently. Namely, the entire top-iy loop is replaced by an effective 
operator insertion. Technically this amounts to performing the integral over its obvious 
internal momentum, after having expanded top and W propagators in remaining momenta 
and M\Y. We note that contributions of internal momentum in the t — W loop of the order 
much less than rrit are suppressed; this is obvious in the unitary gauge which we use in this 
calculation. 

Contributions of diagrams containing the top quark are finite and independent of mt 
for both up and down type quarks (throughout this paper we adopt the approximation in 
which we neglect terms suppressed by inverse powers of quark masses). In particular, we 
have found no terms enhanced by mf, like m^/M^. 

The question which now arises is whether the effective Fermi theory is sufficient for the 
calculation of the remaining diagrams, without the top quark. This turns out not to be 
the case. In diagrams with an external down quark we have to include the contributions 
depicted in Fig. ^ with the top quark replaced by up and charm quarks. Similarly, for an 
external up quark, diagrams like in Fig. ^ contribute, with t — > c. In these diagrams the 
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internal momenta of the order of give contributions to the EDM which are no longer 
cancelled by the corresponding top diagrams. They have to be added to the (divergent) 
effective Fermi theory results. In the sum divergences cancel and logarithms of are 
combined with those of the mass of the second heaviest quark, m^. In this way we arrive at 
our final formulas for the down quark 



e 



and for the up quark: 
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where 0{m? / M'^) denotes terms with an additional suppression by quark or W boson masses; 
Lab = liil'^a/'^^b); Cs is the Riemann zeta function ((^3 = 1.202 . . .), and 6 denotes the CP 
violating invariant 



6 = SIS2S3C1C2C36 



(7) 



where for the CKM matrix we use the convention of p2| , ^ . From our expressions it is not 
obvious that the quark EDM vanishes if any two up- or down-type quarks have degenerate 
masses. This is because we have assumed a hierarchy of masses, and retained only the 
leading terms in the expansions in mass ratios. Such expansion, while delivering a simple 
and compact formula, obscures the symmetry of the full result. 

It is instructive to compare our findings with the results obtained in the Fermi theory 
p^ . In that reference the terms with three powers of logarithms were calculated and if the 
top quark mass is replaced by mass of the W boson they coincide with our Ll^L^b (for dd) 
and —Ll^Lwb — LbcLlg (for du)- However, it turns out that the terms with the highest power 
of logarithms are not the dominant components of the complete result. In fact, the large 
negative coefficient of the single logarithms Lbc leads to more than a complete cancellation of 
the triple log contribution to dd- Partial cancellation is also found in d^- For the numerical 
estimates we use the following values of parameters: 5 = 5 x 10~^, = 0.2, = 0.2 GeV, 
nic = 1.5 GeV, = 4.5 GeV, Mw = 80 GeV. We find 



d. = -0.7 X 10 



_,2 md 



GeV 



e cm 



d^ = -0.3 X IQ-^'^—^ e cm 
GeV 



(8) 



Previous estimates were based only on the terms with three powers of logs; from those we 
find 
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(id(triple log approx.) = + 5 x 10~^^ ^ ^ e cm 

(i^j (triple log approx.) = —0.4 x 10~'^^ ^ ^ e cm (9) 

We see that the "non-leading-log" terms suppress the quark EDM, especially strongly for 
the d quark where they even change the sign of the effect. For the experimentally interesting 
neutron EDM, the d quark EDM is the more important quantity. The suppression factor we 
found decreases the part of the standard model contribution generated by the quark EDM 
and renders it definitely unobservable for the experiments in the near future. 

Inserting current quark masses for rriu^d in eq.(|^) we obtain the following numerical values 

da = -0.7 X 10"^^ e cm for = 10 MeV 

du = -0.15 X 10"^^ e cm for m„ = 5 MeV (10) 

We should stress here that we have neglected terms with more than one power of light 
quark masses m„ rf. They tend to result in an additional suppression, especially for du, where 
ml should be replaced by — rri^ in most terms. 

Regarding the accuracy of our result, it is clear that the issue of what light quark masses 
should be used is most important and requires further study. This is connected with the 
issue of the scale of the running as which is difficult to address since our calculation is of 



the leading order in QCD. For the triple- logarithmic term it is known from p2[ that the 
characteristic gluon virtuality is between c and b quark masses (for dd)- For the non- leading 
terms this scale could be different. Also, since dd results from large cancellations among 
terms which are 5 to 10 times larger than the final value, a precise prediction would require 
taking into account corrections of the order M^/rri^ ^ 20% and ml/ml ^ 10% which can 
change the result by a factor of 2 or so. It is, however, rather certain that the non-leading- 
log terms we presented in this paper significantly reduce the value of dd compared with 
previous estimates. If the next generation of experiments detects non-zero dn it will have to 
be attributed to a different mechanism than the standard model light quark EDM. 
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FIG. 1. Flavor structure of diagrams contributing to quark EDM. Gluon lines are not indicated. 
Here and in the following figures it is understood that an external electric field can interact with 
any charged particle in the diagram. 




(a) (b) 
FIG. 2. Effective Fermi theory diagrams in a model with all quarks lighter than the W boson. 
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FIG. 3. Heavy top quark contributions to da- 
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